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Abstract. We studied the effect of disorder on the superconducting properties of 
polycrystalline MgB2 by specific-heat measurements. In the pristine state, these 
measurements give a bulk confirmation of the presence of two superconducting gaps 
with 2Ao/A;bT' c = 1.3 and 3.9 with nearly equal weights. The scattering introduced by 
irradiation suppresses T c and tends to average the two gaps although less than predicted 
by theory. We also found that by a suitable irradiation process by fast neutrons, a 
substantial bulk increase of dH C 2/dT at T c can be obtained without sacrificing more 
than a few degrees in T c . The upper critical field of the sample after irradiation exceeds 
28 T at T^ 0. 
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1. Introduction 

The recently discovered superconductor MgB 2 has attracted a lot of attention, because 
of its unexpectedly high critical temperature, T c , for a phonon- mediated pairing 
mechanism [U El IH 13 IE]- It has been proposed that such a high T c is due 
to the existence of two superconducting gaps |2J Ej, with this claim supported by 
experiments El- It was suggested that the two gaps open on different parts of the 
Fermi surface [21 Oj- One part is three-dimensional, and arises from the it bonding and 
antibonding orbitals; it gives rise to the smaller gap with 2A 7I .q/A;bT c ~ 1.1 — 1.5. The 
second part consists of nearly cylindrical sheets, which arise from the cx-band, and gives 
rise to the larger gap with 2A o . //c B T , c ~ 3.6 — 4.5. 

The anisotropy of the superconducting gaps and that of the Fermi surface should 
lead to an anisotropic upper critical field, H c2 , and this has indeed been reported in 
MgB 2 jTJIH]. The anisotropy is temperature dependent, the ratio of H c2 in the boron 
plane to that along the c-axis being about 3 near T c and about 7 at low temperature. 

An empirical two-band model based on the presence of two different gaps was 
successfully used to fit specific-heat data obtained on different MgB2 samples in zero 
field |9j. In this model, each gap is assumed to follow the usual temperature dependence 
of the BCS theory, whereas the gap ratio 2A //cbT c is allowed to differ from the BCS 
value 3.5. Both gaps are expected to be sensitive to impurity scattering [2]. An 
interesting prediction is that non-magnetic interband scattering will decrease T c in 
the case of two coupled gaps, whereas Anderson's theorem would rather predict that 
scattering is irrelevant to first order for a single gap JU] ■ I n the limit of strong scattering, 
it has been predicted that both gaps will be averaged up to the point where the single gap 
BCS limit is recovered. This should occur for T c = 27 K [2j. This motivated the present 
study where the behavior of both gaps is followed by bulk specific-heat experiments 
while disorder is introduced in the material. We report heat capacity measurements on 
a polycrystalline sample before and after irradiation by fast neutrons. We found some 
suppression of the larger gap after irradiation, whereas the smaller gap remains quite 
robust. 

Finally, the normal-state resistivity and the irreversibility field of MgB 2 were 
reported to be enhanced by proton irradiation [H]. It has also been shown that 
neutron irradiation can increase the upper critical field of this material [T2]. Here 
we show that the upper critical field of our sample increases from about 18 to 28 T 
after irradiation. Compared to the results obtained on thin films by oxygen post- 
annealing [12] or fast quenching jHj techniques, the specificity of this study lies in the 
use of a bulk determination, in addition to magnetotransport measurements performed 
up to 28 T in a Bitter- type magnet. 
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2. Experimental details 

The polycrystalline sample (HP14) was synthesized from Mg (99.8%) and B (99.7%) 
powders at 900°C for 1 hour in a cubic press at 3 GPa. Fig. Ufa) shows the diamagnetic 
transition of the pristine sample as detected by ac susceptibility at 8 kHz and 0.1 G. The 
midpoint of the transition is T c = 37 K. Neutron irradiation was subsequently performed 
at the Triga research reactor in Vienna. The sample was irradiated to a fast neutron 
fluence (E>0.1 MeV) of 10 22 m -2 , and in a second irradiation step a fluence of 2 x 10 22 
m -2 was added. Neutrons induce defects in MgB 2 mainly by neutron capture of 10 B 
followed by the emission of an alpha particle ^2]- Due to the very small penetration 
range of low energy neutrons in MgB 2 , thermal neutrons would induce damage only 
at the surface of the sample, resulting in a highly inhomogeneous defect structure. In 
order to avoid this problem, a cadmium shield was used, which absorbs nearly all low 
energy neutrons. The introduced disorder can be estimated to 1.7 x 10~ 2 and 5 x 10 -2 
dpa (displacements per atom) after the first and the second irradiation, respectively. 

The transition temperature was measured after each irradiation by ac susceptibility 
as shown in Fig. ^ After irradiation, T c shifted to lower temperature, 34.6 K after the 
first and 30.2 K after the second irradiation, and the transition broadened, but remained 
reasonably sharp compared to the result of similar studies jllj . Therefore the disorder 
introduced by neutron irradiation is fairly homogeneous. 

The specific heat was measured before and after each irradiation in the temperature 
range 2 - 50 K by two different techniques. A relaxation calorimeter was used from 2 to 
16 K |6i, and an adiabatic technique from 16 to 50 K |15j . The resistivity of the sample 
was measured by the four-probe method. Data at fields below 17 T were obtained in 
a superconducting magnet by using a dc current reversal technique, with a current of 
2 mA. Additional measurements from 12 to 28 T were performed in a Bitter coil at the 
GHMFL in Grenoble using an ac synchronous detection technique. These measurements 
were performed with a current of 2 mA at 11.7 Hz. 

3. Results and discussion 

3.1. The increase of H c2 

The specific heat of the sample was measured before and after each irradiation in 
different magnetic fields. Fig. El shows the specific heat difference [C(B) — C(14T)]/T. 
The baseline C(14T) represents the normal state specific heat in the temperature range 
shown in Fig. it was smoothed before subtraction. Since the lattice contribution 
cancels after subtraction, Fig. in fact shows the difference between the electronic 
specific heat in the superconducting and the normal state. The amplitude of the specific- 
heat jump in zero field decreases after each irradiation. Its position is shifted in bulk 
by the magnetic field to a lower temperature. The broadening of the transition in 
increasing fields results from angular averaging since our sample is polycrystalline and 
anisotropic. In MgB 2 , the anisotropy of the H c2 is estimated as 2 - 3 close to T c [HI III IS]- 
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Because of the transition broadening due to anisotropy, the onset of the specific-heat 
anomaly (see arrows in Fig. |2J) was used to define T c . This would correspond to the 
transition temperature of a single crystal for H parallel to the boron planes. At H = 0, 
the critical temperature is 37.4 K, 35.8 K and 31.2 K for the sample before, after 
the first and after the second irradiation respectively. Fig. 01 shows the temperature 
dependence of H c2 obtained from both specific-heat and resistance measurements (before 
and after irradiation). Note that the determinations based on the specific- heat jump 
are no longer significant at high fields owing to the smearing of the jump. Nevertheless, 
they unambiguously demonstrate that the average initial slope (dH c2 /dT) Tc increases 
after each irradiation. 

A positive curvature (PC) of H c2 near T c can been seen in Fig. El at different stages 
of irradiation. This can be attributed to the existence of two superconducting gaps in 
the sample, which may be considered as a particular case of gap anisotropy jTBl ITTj . 
It is pointed out by theory (TTJ HH] that the PC disappears in the "dirty" limit, where 
very strong interband scattering suppresses the anisotropy. In order to investigate the 
change of PC after irradiation, we used the following formula to fit our data near T c up 
to 4 - 5 T: H c2 {T) = H* 2 (l — T/T c ) a [IE], where T c is the transition temperature at zero 
field, H* 2 and a are fitting parameter. We found a within the range of 1.4 ± 0.1 before 
and after each irradiation, and no significant trend was observed. This is consistent 
with the observation that both gaps can still be identified after irradiation, confirming 
that the interband scattering caused by disorder in this sample is not strong enough to 
reach the "dirty" limit, as discussed in Section 3.2. 

In order to verify the increase of H c2 after irradiation, we performed resistance 
measurements in magnetic field up to 17 T on the pristine sample and up to 28 T after 
the second irradiation. The insets of Fig. 01 show typical field and temperature sweeps 
after the second irradiation. For the pristine sample, H c2 (0) is estimated to be about 
18 T, in agreement with the results obtained on single crystals for the field parallel to 
the boron plane After the second irradiation, the transition becomes broader at 

low temperature, as it was already the case for the pristine sample, presumably again 
as a consequence of anisotropy |7J. At T = 1.5 K, the transition is not complete even 
at 28 T, indicating H c2 (0) higher than 28 T. Both specific-heat onsets and resistance 
transition midpoints are plotted in the H — T phase diagram of Fig. 01 The results of 
both determinations of H c2 (T) agree reasonably well in the low field range, the onset 
of calorimetric transitions appearing as an upper limit as expected. The quasi-linear 
increase of H c2 (T) down to T = differs remarkably from the horizontal slope at T — > 
found for conventional superconductors. An explanation for this unusual behavior was 
provided within the two band model JH]- Note that a similar PC was found in several 
borocabides |18j . not necessarily implying the same physical origin |19j . 

For a type-II superconductor in the dirty limit, H c2 (0) oc 7 n pT'c, where 7„ is the 
Sommerfeld coefficient in the normal state, and p is the normal-state resistivity at low 
temperature [20! • After irradiation, T c decreases by 20%, so that the main parameters 
to be investigated for changes are 7„ and p. As will be shown later, 7„ does not change 
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significantly. 

The residual resistivity depends sensitively on disorder. Before irradiation, just 
above T c , p(45K) = 4.0 //O.cm and the residual resistance ratio is RRR = 
p(300K)/p(45K) = 3.3. To estimate the electronic mean free path, we use the following 
formula: 

/=[^] 2 - — , (1) 

e 7 n f F 

where a is the normal-state conductivity and Vp is the Fermi velocity. With 7„ = 
0.14 mJ/K 2 cm 3 [21] and v F = 4.8 x 10 7 cm/s [22], this yields I w 270 A. Since 
/ ^ (£) ~ 42 A [2B|, the sample is initially rather in the clean limit. After the second 
irradiation, the low-temperature resistivity increases to p(45K) = 22.6/xfi.cm (RRR = 
1.4), which gives I ~ 50 A. The sample is now in the intermediate case between the 
clean and dirty limit. In the dirty limit, the upper critical field would be determined 
by the geometrical average between the coherence length and the mean free path 
£ = v^o^- Therefore, at least qualitatively, one is led to conclude that disorder induced 
by irradiation affects H c2 through the reduction of the mean free path. However, we 
note that the sample approaches the dirty limit only after the irradiation, so that the 
discussion given here remains essentially qualitative. Moreover, the two-band nature of 
MgB2 complicates the determination of the mean free path [2*o] . 

Fig. 0] shows the specific heat below 16 K in different magnetic fields at different 
stages of irradiation. At the highest field accessible to our measurements, fioH = 14 
T, C /T reflects nearly the normal-state specific heat. At first view, this would not 
seem to be the case, since fioH c2 (0) is 18 and 28 T before and after the second 
irradiation, respectively. However, it was shown by several experiments that 7 versus 
H already saturates at about H c2 /2 in polycrystals [HI El], so that no significant change 
is expected above 14 T. This is also illustrated by the behavior of the pristine sample 
(see Fig. Ufa)), where specific heat saturates near 8 T. In the inset of Fig. HJ we show 
that the low temperature specific heat at 14 T remains unchanged at different stages of 
the irradiation, indicating that the normal state Sommerfeld coefficient does not change 
(the downturn that occurs at very low temperature in the C/T versus T 2 plot in the 
inset of Fig. HJis believed to be due to the presence of residual magnetic impurities). 
Therefore, as already mentioned above, the variation of H c2 is not due to a change in 
7„. We thus conclude that irradiation merely enhances the scattering, causing both a 
decrease of T c and an increase of the upper critical field. 

These results suggest that other methods able to introduce disorder, such as fast 
quenching or chemical doping, could also serve to increase H c2 , as was realized in A15 
compounds [25]. Indeed, a large increase of H c2 was obtained by metallurgical heat 
treatment of MgB2 thin films fBl E|- The present study shows that this increase of 
H c2 is a bulk property, at least at low fields, where the specific- heat jump is clearly 
observable. It is confirmed by direct resistance measurements at low temperatures and 
high magnetic fields, rather than by extrapolation beyond 10 T in previous work. 
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3.2. Effect of irradiation on the gaps 

As discussed above, irradiation does not reduce the density of states, but merely 
increases the scattering. For a one-band, isotropic superconductor, scattering due to 
non-magnetic impurities would account for the increase of the normal-state resistivity, 
but would not affect T c to first order. For a two-band superconductor, intraband 
scattering would preserve the two-gap structure even for a high concentration of 
impurities [26] . In the present case, a clear suppression of T c is observed; it can only be 
explained by interband scattering between the a- and n- bands EH]- Therefore, 
a distinction should be made between samples for which impurities affect only the 
normal-state resistivity, leaving T c unchanged j2E], and the present irradiated samples, 
for which both T c and resistivity are affected. These distinct situations reflect different 
balances between intra- and interband scattering. If we believe that strong interband 
scattering exists for the present samples, which is rather unusual for MgB 2 [26J, we 
expect changes in the two-gap structure after strong irradiation. We proceed to study 
this point experimentally. 

The electronic specific heat C c /T versus T at zero field is plotted in Fig. 
The lattice contribution determined at 14 T was smoothed and subtracted. The low 
temperature excess of the specific heat with respect to a one-gap BCS curve (hatched 
area in Fig. EJ, which results from the existence of the smaller gap [Oj, becomes less 
pronounced after irradiation. A satisfactory fit can be obtained for the unirradiated 
sample using a two-gap model E]- The fit gives Ao = 6.1 meV for the larger gap 
and 2.1 meV for the smaller gap, with nearly equal weights. Similar fits can be made 
after irradiation. The gap values at each stage of irradiation are given in Table ^ After 
the second irradiation, the quality of the fit was degraded by the broadening of the 
transition. A satisfactory fit was obtained by shifting T c slightly above the midpoint of 
the jump and by lowering 7 n by 2.5%. 

An obvious effect of the irradiation is the suppression of the larger gap, which, 
in a two-gap model, is reflected in the specific heat jump at T c 9J. This gap, which 
originates from the 2D a-band, seems to be more sensitive to defects. The smaller gap, 
which is reflected in the low temperature hump in the specific heat near 10 K, appears 
to be quite robust in absolute value, remaining within 10% of the average 2.15 meV 
in all cases. However, its reduced value 2A //c B ^c increases after each irradiation, due 
to the decrease of T c . Liu et al. predicted that interband scattering in MgB 2 should 
average both gaps, which would finally merge into a single BCS gap when T c reaches 
~ 27 K 2j. The observed gap values are plotted as a function of T c in the inset of Fig. El 
This plot suggests that the single-gap limit will only be reached for T c smaller than the 
anticipated 27 K, although additional irradiations will be required to prove it. From the 
fits, one can also see that the weight of both gaps remains almost unchanged after the 
irradiation. 

Golubov et al. calculated the effect of interband scattering on specific heat due 
to non-magnetic impurities within the two-band model ^0], and pointed out that 
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the impurity scattering should increase the AC/7 n T c ratio, whereas the specific heat 
jump AC should remain nearly constant. According to our experiment, AC decreases 
after irradiation and AC/7 n T c decreases only slightly. This disagreement may be of 
experimental origin and due to the transition broadening. 

The coefficient of the mixed-state electronic specific- heat, 7(H) , provides 
independent information on both gaps. Its highly non-linear increase is shown by the 
plot C(H,T)/T at T < T c (here, 3K) on a logarithmic field scale (see Fig. [HJ). The 
curves at different stages of irradiation lie nearly parallel to each other below 10 T, 
and saturate at the same value of 7 n at 14 T. The small maximum at 8 T for the non- 
irradiated sample is believed to be an artifact due to a residual Schottky contribution 
of magnetic origin. The fact that the three curves are parallel to each other implies 
a simple relation 7 pr istme(-£0 = 7irradiated(aif), where a ~ 1.7 after the first irradiation 
and a ~ 2 after the second one. The inset of Fig. El illustrates this relation: the C/T 
versus T 2 curves at 0.5 T for the pristine sample, 0.84 T after the first irradiation, and 
1 T after the second irradiation point to the same value of j(H, T = 0). 

According to theory |2j and experiment [3 127], both gaps have nearly equal weights, 
i.e. the partial densities of states N a and N n of the associated bands, and the partial 
normal-state Sommerfeld coefficients 7 nCT and 7 n7r , are nearly equal. Recent scanning 
tunneling microscopy measurements on a single crystal show that already in a magnetic 
field of 0.2 T, there is a significant overlap of vortex cores in the 7r-band vortices [28J. 
Very recent specific- heat measurements on a single crystal have revealed that ^y n (H) 
increases very fast with the field, and saturates above a crossover field [XqH = 0.4T, 
whereas J a (H) saturates at 3.3 T (field along the c-axis) or 18 T (field in the ab- 
plane) [29|. It follows that at low fields(i? < 0.5 T) the contribution from the smaller 
gap dominates the behavior of ^(H), in agreement with theory [SU]. On the contrary, 
at high fields (B > 1 T), the contribution from the smaller gap saturates, so that the 
variation of 'y(H) comes from the larger gap. The crossover near 0.4 T can be considered 
as the virtual upper critical field H c2n for the smaller gap. The parallelism of the curves 
in Fig. IHJat high fields (1 - 8 T) provides the proof for the increase of H C 2 associated with 
the larger gap. This is independently measured by the shift of the specific-heat jump 
and resistive transitions. Alternatively, at low fields (B < 0.5 T), below the saturation 
of the contribution 7^ (-£/") associated with the smaller gap, one may deduce from the 
approximate relationship ^(H) = r y nn H/H C 2 that H c27T has increased after irradiation 
by the same factor a as H&a- This observation tends to support the interpretation 
that the superconductivity in the 7r-band is not independent, but rather induced by 
superconductivity in the a-band (HOI EH] . At this point we must recall that our sample 
is a polycrystal, so that information on the anisotropy is lost in the high field region 
dominated by the large gap of the cr-band. More detailed information on this unusual 
system would require irradiation experiments on single crystals. 
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4. Conclusion 

We studied irradiation effects on a polycrystalline sample of MgB2. Irradiation by 
fast neutrons provides a good tool to introduce mterband scattering. Specific-heat 
measurements show a suppression of the larger gap after irradiation, whereas the smaller 
gap A n remains nearly unchanged. In terms of reduced gap value 2A //cb^c, the 
gaps tend to converge, but slower than predicted by considerations of the interband 
scattering. The two-gap feature remains quite robust. We also find that the upper 
critical or crossover field H c27T associated with the smaller gap follows the behavior 
of that of the larger gap after irradiation. Together with resistivity measurements, 
we show that a substantial increase of H c2 can be obtained by irradiation, without 
sacrificing more than a few degrees in T c . If critical currents do not collapse, the 
potential of MgB2 tuned by disorder is worth further investigations in view of high 
field superconducting applications. 
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T C (K) 
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30.2 


7„ (mJ/K 2 .gat) 


0.83 


0.83 


0.81 


A 7r (meV) 


2.07 ±0.20 


2.20 ±0.20 


2.22 ±0.20 


2A 7T /k B T c 


1.3 


1.5 


1.7 


A CT (meV) 


6.21 ±0.30 


5.30 ±0.25 


4.68 ±0.35 


2A a /k B T c 


3.9 


3.6 


3.6 




0.5:0.5 


0.55:0.45 


0.55:0.45 


AC(mJ/K.gat) 


27.6 ± 1 


25.3 ± 1 


20.3 ± 1 


AC/ ln T c 


0.9 


0.85 


0.85 


Table 


1. Gap parameters of MgB 2 


at different stages 


of irradiation obtained by fi 



of the specific heat. 7„ is the normal-state Sommerfeld coefficient. A w and A CT are the 
gap values for the smaller and larger gap. AC is the specific-heat jump. j n7r and 7 nCT 
are the partial Sommerfeld coefficients associated with the tt- and cr-band, respectively; 

Inn + 7ncr = 7n- 
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Figure 1. Diamagnetic transition of the sample measured by etc susceptibility, real 
(x') and imaginary (x") part: (a) before irradiation; (b) after the first irradiation; (c) 
after the second irradiation. 



Figure 2. Transitions at T c in different magnetic fields measured by specific heat, 
AC = C(B = 0) - C(B = 14T): (a) before irradiation, from right to left B = 0, 0.5, 1, 
2, 3 T, (b) after the first irradiation, from right to left B = 0, 0.84, 1.7, 2.5, 3.4, 4.2, 
5.1 T, (c) after the second irradiation, from right to left B = 0, 0.5, 1, 2, 3, 4, 5, 6, 8 
T. Note that the sample was cut and its mass reduced for the irradiations, explaining 
the increased scatter. 



Figure 3. Upper critical field H C 2(T) before and after irradiation: A and full line, 
before the irradiation, determined by the onset of the specific-heat jump; □, before 
the irradiation, determined by the midpoint of the resistance step, o and full line, after 
the first irradiation, specific heat onset; o and full line, after the second irradiation, 
specific heat onset; •, after the second irradiation, transition midpoint in resistance. 
Inset: resistance measurements after the second irradiation: upper right: magnetic 
field sweeps at constant temperature, from left to right: T = 28, 24, 20, 16, 12, 8, 4, 
1.5 K; lower left: temperature sweeps at constant field, from left to right: B = 14, 
12.5, 10, 8.5 T. 



Figure 4. Total specific heat below 16 K in a C/T versus T 2 plot: (a) before 
irradiation, from bottom to top: B — , 0.05, 0.1, 0.2, 0.3, 0.5, 1, 2, 4, 8, and 
14 T; (b) after the first irradiation, from bottom to top: B = 0, 0.17, 0.5, 0.84, 1.7, 
3.4, 6.8, 14 T; (c) after the second irradiation, from bottom to top: B = 0, 0.2, 0.5, 
0.75, 1, 2, 4, 8 T. Inset: low temperature specific heat at 14 T before and after each 
irradiation (three data sets) . The extrapolation to T = shows that j n does not vary. 



Figure 5. Electronic specific heat C c /T versus T with two-gap fits from top 
to bottom: before irradiation, after the first irradiation, after the second irradiation, 
respectively. The latter two curves are shifted by 0.25 mJ/K 2 gat for clarity. The 
dashed line represents the BCS single-gap model with the same T c and 7„ as the 
sample after the second irradiation. Inset: variation of both gaps as a function of T c 
(note that the scale is reversed to show increasing scattering from left to right). The 
dashed line is the isotropic BCS limit. The star * is a calculated point where the gaps 
should converge to the BCS value according to Ref. 0- 
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Figure 6. Specific heat C/T at 3 K versus magnetic field on a logarithmic scale, 
o before irradiation, • after the first irradiation, A after the second irradiation. Inset: 
from bottom to top, low temperature specific heat at 0.5 T before irradiation, 0.84 T 
after the first irradiation and 1 T after the second irradiation, respectively. 
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